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Abstract 

Multiple zeta values (MZVs, also called Euler sums or multiple har- 
monic series) are nested generalizations of the classical Riemann zeta func- 
tion evaluated at integer values. The fact that an integral representation 
of MZVs obeys a shuffle product rule allows the possibility of a combi- 
natorial approach to them. Using this approach we prove a longstanding 
conjecture of Don Zagier about MZVs with certain repeated arguments. 
We also prove a similar cyclic sum identity. Finally, we present extensive 
computational evidence supporting an infinite family of conjectured MZV 
identities that simultaneously generalize the Zagier identity. 

1 Introduction 

In this paper, we continue our study of multiple zeta values (MZVs), sometimes 
also called Euler sums or Zagier sums, defined by 

fe 

C(si,...,s fc ) := U n j Sj 

ni>n2>..->n k >0 j=l 

with Sj 6 Z + and si > 1 to ensure the convergence. The integer k is called the 
depth of the sum C(si , . . . , Sfc). 

MZVs can be generalized in many ways. In particular, they are instances of 
multidimensional polylogarithms [|| . Such sums have recently attracted much at- 
tention, in part since there are many fascinating identities among them. The ap- 
plications of MZVs involve some unexpected fields, such as high energy physics 
and knot theory — see for a list of references. 

Hoffman in his study |6) of the ^-product of MZVs (which we call the "stuffle" 
product in |2|) distinguishes between "algebraic" and "non-algebraic" relations 
among MZVs — the latter ones involve a limiting process in some essential way. 
In the same spirit we note that some non-trivial MZV identities are conse- 
quences of discrete (combinatorial) relationships involving the shuffle product. 
Hints that this may be the case include the occurrence of binomial coefficients 
(e.g., in (|l0|)). In the present paper we follow the combinatorial approach by ex- 
ploring the combinatorial content of the shuffle product rule (^|) for the integral 
representation [pi of MZVs. 

In Section |2jwe list some factorial identities on which we base our later re- 
sults. In Section || we introduce the shuffle algebra and in Section || we prove 
some combinatorial identities holding in this algebra. The relevance of the shuf- 
fle algebra for studying MZVs originates in the iterated integral representation 
of MZVs which we briefly recall in Section |[ In Section |we use shuffle iden- 
tities to prove the longstanding conjecture of Don Zagier [0 [j], |[: 
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(where the notation {A}™ indicates n successive instances of the integer se- 
quence X), as well as the similar "dressed with 2" identity: 



(4n + 3)!' 



where s runs over all 2n + 1 possible insertions of the number 2 in the string 
{3, 1}". Finally, in Section @ we present extensive numerical evidence for our 
new conjecture, which in a rotationally symmetric way generalizes (by insertions 
of groups of 2's) the Zagier identity. For an illustration, one very simple instance 
of our conjecture reads 

7T 10 

C(3, 2, 2, 1, 2) + C(2, 2, 3, 2, 1) + C(2, 3, 1, 2, 2) = — . 

2 Factorial Identities 

In the main part of the paper we will require the following identities. The 
proofs are easy by any of several methods (generating functions, WZ theory, 
etc.); therefore we skip them. 

Lemma 1 For any non-negative integer n we have 

(-iy 2 2n+1 



E (2n + 2r+l)!(2n-2r+l)! (4n + 2)!' 



Lemma 2 For any non-negative integer n we have 

<A (-l) r (2r + l) 4 n 



(2n + 1 - 2r)!(2n + 3 + 2r)! (4n + 3)!' 



Lemma 3 For any non-negative integer n we have 

STf l\rfo„ , i\f 2n+1 \ _ / 1 tfn = 



r=0 



(1) 



3 The Shuffle Algebra 

Let A denote a finite alphabet (set of letters). By a word on the alphabet A 
we mean a (possibly empty) sequence of letters from A. By A* we denote the 
set of all words on the alphabet A. For w G A*, let w k denote the sequence of 
k consecutive occurrences of w. A polynomial on A over Q is a rational linear 
combination of words on A. The set of all such polynomials is denoted by Q(^4). 
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On Q(A) we introduce the binary operation lu ("shuffle product"), which 
is defined, for any u,v £ A* (u — x\ . . . x n and v = x n +i . . . x n + mi Xk G A for 
1 < k < n + m) by 

uujv := y^a: .(i)a; .(2) ■ ■ .x„( n+m ), (3) 

where the sum is over all ( n ~^ m ) permutations a of the set {1, 2, . . . , n + m} 
which satisfy c^ 1 ^') < cr _1 (fc) for all 1 < j < k < n and n+l < j < k < n + m. 
In other words, the sum is over all words (counting multiplicity) of length n + m 
in which the relative orders of the letters x\, . . . ,x n and £n+ij . . . , x n+m are 
preserved. The definition (||) extends linearly on the entire domain Q(^l) xQ(^l). 

Example. Let A = {A, B}. In Q{A) we have 
2AB lu (3BA - AB) = 12AB 2 A + 12BA 2 B + 2(AB) 2 + 6(BA) 2 - 8A 2 B 2 . 



4 Identities Involving Shuffles {ABf with (AB) q 

Throughout the rest of the paper we assume that the alphabet A contains 
exactly two letters A and B. 

Definition 1 Let p, q and j be non-negative integers subject to min(p, q) > 
j. Let Sp+qj denote the set of those words occurring in (AB) P lu (AB) q that 
contain the subword A 2 exactly j times. 

Definition ^ is sound, since the set S p + q .j is the same for any partition of the 
number p + q into two parts as long as both parts are greater than or equal to j. 
This would of course not be true if we instead considered the full expansion of 
(AB) P lu (AB) q (that is, counting the multiplicity of words): see Proposition [t] 
in which we calculate these multiplicities explicitly. 

Side remark. The set S p + q j has cardinality ( P 2j 9 ) ■ Indeed, any word in 
Sp+qj can be considered to be partitioned into p + q consecutive blocks of 
length 2. Clearly, the locations of the subwords A 2 and B 2 are consistent with 
this partitioning. Since there are j blocks containing A 2 , they must be interlaced 
with another j blocks containing B 2 , and the choice of the positions of these 
j + j = 2j blocks together with the shuffle rule W) determines the rest of the 
word in question. Therefore there are exactly ( p ^y) elements in S p + q .j. 

Definition 2 Let p, q, j be as in Definition [J. By T p+q j we will denote the sum 
of all words in S p + q j . 



Proposition 1 For any non-negative integers p and q we have 

min(p,g) . _ O \ 

{ABfm{ABY= £ V ■("■ '' .-' • /;,,,, 

j=o 
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Proof. Let u be an arbitrary but fixed word from S p + q j . Let us see how many 
times u arises in (AB) P lu (AB) q . This is the same as counting in how many 
ways the letters of u can be colored in two colors (blue letters coming from 
(AB) P and red letters coming from (AB) q ) in a coloring that is consistent with 
the shuffle rule @. 

There are p + q A's in u, of which 2j A's are contained in factors A 2 and 
p + q — 2j A's are surrounded by B's from both sides (or possibly from one 
side if we are looking at the leading A). Of the latter p + q — 2j "single" A's, 
p — j are colored blue. Thus the coloring of the single A's contributes a factor 
of ( p +t} 2i ) to the multiplicity of u in (AB)p lu (AB)«. There are exactly j 
factors A 2 (and thus exactly j factors B 2 ) in u, each of which can be colored 
in two ways (blue-red or red-blue), thus contributing a factor of 2 J • 2 J = 4 J to 
the multiplicity of u in (AB) P lu (AB) q . What remains to do is to color the 
"single" B's, whose coloring is now determined uniquely by the choices made so 
far, together with the shuffle rule (||). □ 

Corollary 1 For any non-negative integer n we have 

n 

(-l) r [(AB) n -' r lu [AB) n+r ] = A n (A 2 B 2 ) n . (4) 

r — — n 

Proof. Using Proposition the left-hand side of (0) is equal to 

n min(n-r,n+r) . . 

^— ' / r~i \n — r — 1 1 

r=— n j=0 x 7 

which after reordering is 

j=0 r=j— n ^ ' 

Putting N :— n — j in the inner sum turns it into 




which is equal to 1 if N — (i.e. j = n) whereas for N > (i.e. j < n) it is a 
disguise of (1 - 1) 2N ■ (-l) N which is 0. Thus, (|) is equal to 4" • T 2n>n which 
is indeed the right-hand side of ([|), and the proof is finished. □ 
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Corollary 2 For any non-negative integer n we have 

n 

^(-l) r (2r + l) [{AB) n - r lu [AB) n+1+r ] = 



r=0 



4" • ^{A 2 B 2 ) r AB{A 2 B 2 ) n - r + Y^{A 2 B 2 ) r - l A 2 BAB 2 (A 2 B 2 ) n - r . (6) 

\r=0 r=l / 

Proof. First we show that, for any non-negative integer n, we have 

n 

^(-l) r (2r + 1) [(AB) n - r lu (AB) n+1+r ] = 4"T 2n+1 , n . (7) 

r=0 

As in the proof of Corollary [l] we proceed in three steps: (i) evaluating the 
shuffle products by Proposition [I], (ii) swapping the sums, (iii) doing the inner 
sum. 

Using Proposition GJ, the left-hand side of (fjj) can be written as 

r=0 j=0 ^ •' ' 

which after reordering is equal to 

which by Lemma || (with n — j in the place of n) is equal to 4™T 2n+ i ira . 

Now T 2n+ i i „ is the sum of words arising in the shuffle (AB) n lu (AB) n+1 
and containing n factors A 2 and n factors B 2 . Thus, there is exactly one single 
A and exactly one single B, which clearly have to be adjacent, and thus forming 
a factor AB or BA. In the parentheses on the right-hand side of (||), the first 
summand accounts for those summands from T2 n +i. n that contain AB, while 
the second summand accounts for those summands from T2„+i. n that contain 
BA. This completes the proof of (||) . □ 



5 Integral Representation of MZVs 

Let us recall that we are working with the alphabet A = {A, B}. Throughout 
the rest of this paper we identify the letter A with the differential form dx/x 
and the letter B with the differential form dx/(l — x). 

The MZV C( s i: • ■ • , Sfc) admits the (s\ + s 2 H — • + sjj-dimensional iterated 
integral representation 

C( Sl ,...,s k )= [ A*- 1 BA»- 1 B--.A'*- 1 B, si>l. (8) 
Jo 
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The explicit observation that MZVs are values of iterated integrals is apparently 
due to Maxim Kontsevich H]. Less formally, such representations go as far 
back as Euler. The representation (||) is a very special instance of the iterated 
integral representation of multidimensional poly logarithms @J — see there for the 
exact definition of the iterated integral (|J) , which however is not critical for our 
purposes. 

Indeed, the only property of iterated integrals that we use in this paper is 
that their products obey the "shuffle rule," that is |l^, ^ 

' u )il! v ) = Ij u ^ v) (9) 

if we view the products of differential 1-forms in U and V as words in the shuffle 
algebra (Section ||). Clearly, (||) motivated our interest in shuffle identities 
(Section |). 

An intriguing aspect of (0) is the bridge between analytical (transcendental) 
and discrete nature of MZVs. Although the present paper deals only with MZVs, 
the ideas used here are applicable to more general nested sums (alternating sums, 
multidimensional poly logarithms @1) since, as already mentioned above, these 
sums admit integral representations which generalize (|^). 

Example. We provide a combinatorial derivation of Eulcr's decomposition 
formula (s, t > 2) 



3=1 V S ' 



t 

E 



s + t-j 

t-j 



<Z(s + t-j,j). (10) 



Let us consider the product P := A S ~ 1 B lu A f ~ 1 B. Clearly, any term in P 
must end with a B. The terms in P in which the trailing B comes from the 
A S ~ 1 B operand are accounted for by 



s+t-i _ 1 

l K 1 \ Ak-l o As+t-k-l 



^ \t-l 

k=t 



A* -1 BA a+ *-* -1 B, (11) 



with the binomial coefficient counting the number of ways in which all A's from 
the A t ~ 1 B operand can be inserted in the leading block of A's in the shuffled 
string. Similarly, those terms in P in which the trailing B comes from the 
A t ^ 1 B operand are accounted for by 



s+t-l , k _ 1 



k=s 



E [ i)A k - 1 BA'+ t - k - 1 B. (12) 



Summing up ( |ll| ) and (|12|), substituting k := s + t — j and using (|9||8j) gives 

Tol). 
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6 Proof of the Zagier Conjecture 

From Section [l] we recall that, in the context of integer sequences, we use the 
notation {X}™ to indicate n > successive instances of the sequence X. 

Theorem 1 (The Zagier Conjecture) For any positive integer n we have 

2-7T 4 ™ 

«^ 1 > W )=(^" (13) 

Proof. Using (^,||) , Corollary [j] implies 

n 

£ (-l) r C({2}"- r )C({2}" +r ) = 4"C({3,1}"). 

r— — n 

Application of the evaluation 

C({2D = (U) 



which was proven in U , gives 

4"C({3,l} n )^ 4 " £ 



(-iy 



^ (2n-2r + l)!(2n + 2r+l)! 

r=—n v ' x ' 

which by Lemma [l] is equivalent to 

4"C({3,1}") = 7T 4 " 



(4n + 2)! 

After dividing the last equation by 4™ we get (fL3|). □ 



The first proof of (|l3|) appears in ||. It may be viewed as the first non- 
commutative extension of Euler's evaluation of C(2n). 

Theorem 2 Let n be a "positive integer, and let I denote the set of all 2n + 1 
possible insertions of the number 2 in the string {3, 1}". Then 

_-4ri+2 

z«">-wm- <15) 

sEl K ' 

Proof. Using (|9|,||), Corollary || implies 

n 

£(-1)^ + l)C({2}»-'-)C({2}" +1+ '-) = 4" £ C(S). 
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Indeed, the first term in the parentheses on the right-hand side of (|^) translates 
to J2r=o C({3, l} r , 2, {3, while the second term translates to 

Er=iC({3,l} r - 1 ,3,2,l,{3,l}"^)- Application of <M) gives 



(-l) r (2r + l) 

-> (2n - : 

which by Lemma || is equivalent to 



4™ V C(s) = 7T 4 "+ 2 V — 

1 ^ (2n-2r+l)!(2n + 3 + 2r)! 



4"^C(s)=^ 4 " +2 — 



(4n + 3)! 

sei v ' 

After dividing the last equation by 4™ we get ([L5j). □ 



7 Conjectured Generalizations of the Zagier 
Identity 

To notationally ease our generalization, we define 

Z(m , . . . , m 2 „) := C({2}™°, 3, {2} m \ 1, {2}™ 2 , . . . , 3, {2} m --\ 1, {2}™-) , 

(16) 

with {2} m; * inserted after the j-th element of the string {3, 1}". For example, 
Z(2, 0,1) = ((2,2,3,1, 2). 

Conjecture 1 For any sequence S = (mo, . . . ,m 2n ) of 2n + 1 non-negative 
integers, we have 

2n 4n+2M 

!>< tfS >= (4„ + 2M + l)r < 17 > 

3=0 

where M :— X) 2 "o TO * an ^ ^ * s ^ e C U C ^ C permutation operator, that is, 
C j (mo, . . . ,m 2 n) ■= {m 2n -j+i, ■ ■ ■ ,m 2n ,m , . . . ,m 2n -j)- 

Remark. Taking into account ([TJ) we see that the right-hand side of (^) is 
equal to C({2} 2n+M ). 

In Section |^ we proved ([Tt]) for the cases M — and M = 1. For n = 0, 
trivially reduces to the known evaluation (|l4|). If all m^s are equal, ( p~7| ) 
specializes to conjecture (18) of jjj. 

Since MZV duality f| |] implies that 

Z(S) = Z{S), (18) 

where S := (m 2n , . . . , too) is the reverse of S, the conjecture dl^ ) can be also 
reformulated as a sum over all permutations in the dihedral group D 2n+ \. In 
our formulation we sum over the cyclic group C 2n+ \. 
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7.1 Integer Relations 

An integer relation || for a vector of complex numbers z G C" is a non-zero 
vector of integers a € Z n such that 

a\Zi H V a n z n = 0. 

Conjecture [l] was discovered numerically (via its special instances) using the 
PSLQ algorithm for discovering integer relations [[!] and the fast method for 
numerical evaluation of MZVs using the Holder convolution ||. All cases of 
( p~7[ ) with depth 2n + M < 13 were checked numerically at the precision of 2000 
digits. This amounted to checking 747 such identities, even after excluding the 
cases with n = or M < 1, for which proofs have been known before or are 
presented in this paper. 



7.2 Other Conjectured Identities 

For any two fixed integers n, M > 0, let us consider the vector V n> M of values 



Z(mo, mi, . . . , TO2n) defined by (16) and subject to: m, e Z, rm > (0 < 



i < 2n) and Yli=o mi = M. We assume that the entries of V n ,M are listed in 
some arbitrary (but fixed) order, and that of any two Z-terms related by the 
duality (jl^), exactly one is present in V n> M , in order to exclude trivial duplicates. 
Additionally, wc append to V n ,M the value Z{2n + M) := ( ({2} 2n+M ). 

If we restrict our attention to the putative identities of the form (|l7|), then 
the number of (linearly independent) relations of this type can be computed 
via Polya Theory (see, e.g., as the number of orbits in the action of the 
dihedral group £>2n+i on the set of functions / : {0,1,..., 2n} — > N subject 
to 53i=o /(*) = (Let us recall from Section that we have verified ( \n\ ) 
numerically in the range 2n + M < 13.) 

On the other hand, integer relations for ~V n M can be discovered empirically 
using integer relation algorithms, regardless of whether their structure is com- 
patible with (|lj]) or not. In Figure |l| we list, for some modest values of n and 
M, in lightface the number of (|l7j)-type putative relations for V n ,Mi and in 
boldface the number of relations for V n ,M detected empirically using the PSLQ 
algorithm || using the numerical precision of 5000 decimal places. (In both 
cases we count the number of linearly independent relations.) These values (as 
well as some others, not included in Figure |l|) suggest that the scheme ( p7| ) 
exhaustively describes all integer relations for V n .M in the cases when n < 1 or 
M < 2, while in the remaining cases, additional relations were detected. 
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M 


1 


2 


3 


n 








1 


1,1 


2,2 


3, 3 


2 


1,1 


3,3 


5, 6 


3 


1,1 


4,4 


8,10 


4 


1,1 


5,5 


12,15 



Figure 1: Number of linearly independent integer relations for V n ,M- 

Example. Here is a family of identities among Z-values (but not of type 
(^7|)), for which we have extensive numerical evidence: 

Conjecture 2 For any non-negative integers 01,02,03,61,62, we have 

Z(oi, 61, a 2 ,6 2 ,a 3 ) + Z(o 2 , 61, 03,62,01) + Z(a 3 , 61, 01,62,02) 
= Z(ai, 6 2 , a 2 ,6i,o 3 ) + Z(a 2l 6 2 , 03,61,01) + Z(o 3 , 6 2 ,ai,6i,a 2 ). 
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